C<| ' Investigation of next-to-leading effects in CCFM* 

o ■ 

^ ; H. Jung 

Department of Physics, Lund University, Sweden 

(N 

> 

cn 
^ , 

. The CCFM [1] evolution equation in the framework of fci-factorization 

I and its practical realization in the Monte Carlo program Cascade [2] has 

■ been shown to be very successful in describing a bulk of experimental mea- 

, surements [2-5], which were not described in the collinear approach. How- 

Q^' ever, as BFKL, the CCFM equation was derived in the high energy ap- 

proximation keeping only the singular terms (i.e. 1/z and 1/(1 — z)) in the 
splitting function Pg. The question arises, whether the other terms, which 
' are present in the DGLAP splitting function, are already small enough to be 

neglected at the energies of present colliders. Also the scale in the running 
Os was originally treated differently for the small and large z parts. 



The effect of formally next-to-leading contributions to the CCFM evo- 
lution equation are discussed. 



1. Introduction 



X 

. - - . 2. Next-to-leading effects 

The splitting of ki-i — > kiPi, where k (p) are the four-momentum vectors 
of the propagator (emitted) gluon, respectively, with momentum fractions 
Xi-i, Xi and the splitting variable z = Xi/xi-i, is described by the splitting 
function Pg. The original CCFM [1] splitting function Pg was given by: 

1 — z z 

where q = pt/{l — z) and the non-Sudakov form factor A„s was defined as: 

log AUz, q\ kl) = -a, £ ^ / '^Qi.h - q)e{q - z'q) (2.2) 
* Presented at DIS 2002 Hadronic Final States - Cracow 2 May 2002 
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Here only the singular terms 1/z and 1/(1 — z) were included and for sim- 
plicity the scale in the running was not treated in the same manner for 
the small and large z part. 

In the high energy approximation, the inclusion of the non-singular 
terms in the splitting function Pg as well as changes in the scale of ag 
are considered as next-to-leading effects. In the following we investigate the 
numerical importance of these effects at present collider energies. 



2. 1 . Scale of as 

Due to the complicated structure of the CCFM splitting function, for 
simplicity the transverse momentum of the propagator gluon, fc^, was used 
as the scale in the running as, whereas next-to- leading order calculations 
suggest, that the proper scale is the transverse momentum of the emitted 
gluon, pf, for full range of z (for a summary of the arguments see [4]). 






^ : n r d o <r6^^(>^ 








/ q = 10, k = 5GeV 

1 




-2 - 1 
1 1 




z 



Fig. 1. The non-Sudakov form factor A„s for three different values of kt/qt as a 
function of the splitting variable z according to eq.(2.2) (solid) and eq.(2.5) (dotted). 

Changing the scale in ag from kt to pt (with ~ gr for z ^ 0) , also the 
non-Sudakov form factor needs to be properly changed, resulting in: 



(2.3) 



1 — z z 

log A„, = - ['^ [ ^asiq')@{kt - q')@{q' - z'q) (2.4) 
which leads to: 



logA^ 



1 dz' f^l dq 
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z' J(z'qY^ q'"^ \og{q'/kQCD) 



(2.5) 



Due to the angular ordering constraint q' > z'q, q' can become very small 
and even q' < Aqcd at small values of z'. Thus a cutoff is required. For 
z'q < qo = 0.71 GeV we fix 05(^0) = 0.5, but keep the full angular ordering 
constraint in the integral. 
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Fig. 2. The splitting function Pg for three different values ofkt/qt as a function of 
the splitting variable z according to eq.(2.1,2.2) (solid) eq.(2.3,2.5) (dotted), and 
with A„g = 1 (dashed) 

In Fig. 1 we compare the new non-Sudakov form faetor A„s with the 
standard one for three different values of kt/qt- It is interesting to note, 
that everywhere the very small z values are highly suppressed. 

In Fig. 2 the splitting function Pg (dotted) is plotted as a function of z 
for three different values of kt/qt- Also shown for comparison is the split- 
ting function without any suppression from the non-Sudakov form factor 
(A„s = 1 dashed) and the standard version of the splitting function from 
eqs.(2.1,2.2) (solid). One can clearly see how the different non-Sudakov form 
factors suppress the small z region of Pg. 



2. 2. Non-singular terms in Splitting function 

Another source of next-to-leading-log corrections is the gluon splitting 
function itself. At very high energies, the 1/z term in Pgg^ included in BFKL 
and CCFM, will certainly be dominant. However, the question is whether 
including just this term is sufficient at energies available at present colliders. 

The implementation of the full DGLAP splitting function into CCFM 
is problematic. Naively one would simply replace — > — 2 + z{l — z) 
in the CCFM splitting function. But this can lead to negative branching 
probabilities. 

In [4] it was suggested to use: 

P{z, q, k) = a, {kf) (^^^ + (1 - B)z{l - z)^ A„,(2, q, k) (2.6) 

+as((l-^)V) [j^^+Bz{l-z)^ 

where S is a parameter to be chosen arbitrarily between and 1, we take 
B = 0.5. As a consequence of the replacement, the Sudakov form factor 
will change, but also the non-Sudakov form factor needs to be replaced by: 

log = -as (kt) £ dz' + (1 - B)z{l - z)^ j ^Q{k-q')Q{q'-z 

(2.7) 
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Fig. 3. The non-Sudakov form factor Ans for three different values of kt/qt as a 
function of the splitting variable z according to eq.(2.2) (solid) and eq.(2.7) (dotted). 

In Fig. 3 we compare the new non-Sudakov form factor Ans with the 
standard one for throe different vahics of ki / qi . 
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Fig. 4. The splitting function Pg for three different values of kt/qt os a function of 
the splitting variable z according to eq.(2.1,2.2) (solid) eq.(2.1,2.1) (dotted), and 
with Ans = 1 (dashed) 

In Fig. 4 the splitting function Pg (dotted) is plotted as a function of z 
for three different values of kt/qt- Also shown for comparison is the split- 
ting function without non-Sudakov form factor (A„s = 1 dashed) and the 
standard version of the splitting function from eqs.(2.1,2.2) (solid). Here 
the effect of the different form of the splitting function Pg becomes obvious 
already at values of z ~ 0.5, whereas the non-Sudakov form factor is sim- 
ilar to the standard one. One should note that especially in the region of 
medium z, the new branching probability (including the non-singular terms) 
becomes smaller. 

2.3. Consequences for forward jet production 

In Fig. 5 we show the predictions for forward jet production at HERA [6] 
for the different scenarios discussed above. All cases have been re-fitted 
to the structure function F2, with a similarly good 'x^ /ndf. It becomes 
obvious, that the prediction for forward jet production is rather sensitive to 
the details of the gluon splitting function. 
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Fig. 5. The cross section for forward jet production as a function of x, compared 
to HI data [6] 
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